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Abstract 



Infrared and collinear safe event shape distributions and their mean values 
are determined in e"'"e~ collisions at centre-of-mass energies between 45 and 
202 GeV. A phenomenological analysis based on power correction models in- 
cluding hadron mass effects for both differential distributions and mean values 
is presented. Using power corrections, as is extracted from the mean values and 
shapes. In an alternative approach, renormalisation group invariance (RGI) is 
used as an explicit constraint, leading to a consistent description of mean val- 
ues without the need for sizeable power corrections. The QCD /^-function is 
precisely measured using this approach. From the DELPHI data on Thrust, 
including data from low energy experiments, one finds 

(3q = 7.86 ± 0.32 

for the one loop coefficient of the /3-function or, assuming QCD, 

rif = 4.75 ± 0.44 

for the number of active flavours. These values agree well with the QCD expec- 
tation of (3o = 7.67 and nf = 5. A direct measurement of the full logarithmic 
energy slope excludes light gluinos with a mass below 5 GeV. 
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1 Introduction 

The decrease of the strong couphng parameter, cts, with increasing energy, or mo- 
mentum transfer, Q, and the related properties of asymptotic freedom and confinement 
are striking phenomena of Quantum Chromodynamics, QCD, the gauge theory of strong 
interactions. Besides the measurement of the strong couphng itself, the precise measure- 
ment of its energy dependence is an experimental task of fundamental importance. This 
energy dependence is governed by the (5 function, defined as [1] : 



- PM (1) 

The coefficients Pi are given in the Equations 14. 

In principle, the study of event shape observables (e.g. Thrust) in e+e" annihilation 

as a function of energy permits these determinations. Event shape observables, however, 
are obscured by the effects of non-perturbative hadronisation. This influence is expected 
to vanish with increasing energy (going dominantly proportional to 1/E). A similar 
dependence is present in the so-called infrared renormalons appearing in perturbation 
theory [2]. Both phenomena are often considered to originate from the same physics. 

The theoretical analysis of power terms indicates some properties which are univer- 
sal to all observables. A coherent comparison of the power correction models and the 
proposed universality is one topic of this paper. 

Power corrections are subject to ambiguities. They depend on the borderline dividing 
perturbative and non-perturbative physics in the models; this borderline is a matter of 
convention. Usually the perturbative terms are treated in an 0{a'l) approximation, the 
remainder being taken as a power correction. The inclusion of higher order perturbative 
corrections in general will reduce the size of the power terms. 

Power correction models are now considered as established. However, in view of a pre- 
cise measurement of the strong coupling, power corrections and the ambiguity involved 
are obstructive. Quantities for which power corrections are minimised should be empha- 
sised. Moreover the reason for the success of power corrections and their magnitudes is 
"not yet fully understood" [3]. Hence a critical review of other theoretical methods for 
many experimental observables is in order. 

The second focus of the phenomenological analysis presented in this paper is on the 
study of renormalisation group invariant perturbation theory (RGI) [4-6]. Here the prop- 
erty of complete renormalisation group invariance is used, leading to predictions without 
the freedom arising from the choice of renormalisation scheme or scale. The theoretical 
ansatz employed [5], however, only applies to "fully inclusive" observables depending on 
a single energy scale, such as total cross-sections [7,8], or mean values of event shape 
observables. A thorough test of this theoretical method is presented here based on the 
energy dependence of the means of the distributions of seven event shape observables. 
The convincing success of this test implies that the size of the power corrections found 
for the power correction models in the MS renormalisation scheme can be predicted 
using RGI perturbation theory. Furthermore RGI perturbation theory allows a direct 
measurement of the /3-function of QCD avoiding any scheme dependence. 

Since the goal of this analysis is a study of the scale dependence, data with a wide 
range of centre-of-mass energies are needed. Therefore data have been used from the 
high energy running of Lep up to centre-of-mass energies of 202 GeV, from the energies 
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around the Z pole, and also radiative events with a reduced centre-of-mass energy of the 
hadronic system due to hard photon radiation. Additionally for some part of the analysis 
low energy data from other experiments are included. 

The organisation of this paper is as follows: Section 2 discusses briefly the apparatus, 
the data and the analysis of the Lep2 data and the radiative events used to extract data 
below the Z mass scale. The observables used throughout the paper are introduced and 
their dependence on the masses of the final state hadrons is discussed. Finally, proce- 
dures to determine systematic uncertainties are specified. Section 3 presents the data 
on shape distributions and their means and compares them to predictions of prominent 
fragmentation models. Section 4 similarly makes comparisons with analytic power model 
predictions [9-11], together with a comparison of the non-perturbative parameters for the 
different observables. For the first time the prediction of power corrections for the energy- 
energy correlation (EEC) is compared to experimental data and evidence for power shifts 
at the three jet phase space boundary is given. Section 5 contains the power correction 
analysis for mean values of exponentiating event shape observables based on the prescrip- 
tion of [12,13]. Simple power fits [14,15] are then presented for all shape observable means 
and it is shown that the size of the power terms correlates with that of the corresponding 
second order pcrturbative contribution. Section 6 briefly recalls the basics of the RGl 
method as given in [5] and confronts it with the data on shape observable means. It is 
shown that this method describes the data very well. The inclusion of power terms in the 
fit shows no indication of significant non-perturbative effects. Consequently the size of 
the power terms determined in the previous chapter is estimated from the RGI method 
and shown to agree with experimental data. By applying the RGl method to different 
shape observable means as measured by Delphi, especially to the data on (1 — Thrust) 
combined with results of other experiments at low energy, a direct precise measurement 
of the QCD /3-function is obtained. Finally we summarise and conclude in Section 7. 

2 Detector, data and data analysis 

Delphi is a hermetic detector with a solenoidal magnetic field of 1.2T. The tracking 
detectors, situated in front of the electro-magnetic calorimeters are a silicon micro- vertex 
detector (VD), a combined jet/proportional chamber inner detector (ID), a time projec- 
tion chamber (TPC) as the major tracking device, and the streamer tube detector OD in 
the barrel region. The forward region is covered by silicon mini-strip and pixel detectors 
(VFT) and by the drift chamber detectors FCA and FCB. 

The electromagnetic calorimeters are the high density projection chamber HPC in the 
barrel, and the lead-glass calorimeter FEMC in the forward region. Detailed information 
about the design and performance of DELPHI can be found in [16,17]. 

The phenomenological analysis of the event shape data presented in the following 
sections relies on the Delphi data measured at the Z-peak, the off-peak energies of 89 
and 93GeV, as well as the published data between the Z peak and 183GeV[18,19]. In 
addition the data measured at centre-of-mass energies y/s between 189 and 202 GeV and 
from radiative events at mean hadronic centre-of-mass energies of 45, 66 and 78 GeV 
are presented. The number of events accepted in the analysis at these energies and the 
integrated luminosities collected are given in Table 1. 
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- 
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9509. 


17.0 


148 • 10^ 
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14007. 


13757. 


18.6 


230 • 10^ 


0.844 


0.999 


133.2 


291. 


69.2 


11.9 


850 


0.958 


0.999 


161.4 


147. 


32.3 


10.1 




282 


0.711 


0.992 


172.3 


121. 


27.5 


9.9 


219 


0.729 


0.956 




183.1 


100. 


23.4 


54.0 


1035 


0.733 


0.930 


189.2 


99.8 


21.1 


150.7 


2774 


0.749 


0.909 


192.2 


96.0 


20.2 


25.8 




433 


0.757 


0.889 


196.2 


90.0 


19.2 


771 




1288 


0.754 


0.876 


200.1 


85.2 


18.2 


80.9 


1281 


0.767 


0.857 


202.1 


83.3 


17.7 


40.0 


624 


0.764 


0.890 



Table 1: Nominal centre-of-mass energies Ecm, cross-sections, without (a), and with ISR 
cut {(Tcs), luminosities {£), the number of selected events (A^sei), the efficiencies (e) and 
the purities (p). 

2.1 Selection and analysis of high energy data 

In order to select well measured particles, the cuts given in the upper part of Table 2 
have been applied. The cuts in the lower part of the table are used to select e+e" — > 
Z/7 — > qq events and to suppress background processes such as two-photon interactions, 
beam-gas and beam-wall interactions, leptonic final states and, for the Lep2 analysis, 
initial state radiation (ISR) and four-fermion background. 

At energies above 91.2 GeV, the high cross-section of the Z resonance peak raises 
the possibility of hard ISR allowing the creation of a nearly on-shell Z boson. These 
"radiative return events" constitute a large fraction of all hadronic events. The initial 
state photons are typically aligned along the beam direction and are rarely identified 
inside the detector. In order to evaluate the effective hadronic centre-of-mass energy 
a/s^ of an event, taking ISR into account, the procedure described in [20] was applied. 
It is based on a fit imposing four-momentum conservation to measured jet four-momenta 
(including estimates of their uncertainties) . Several assumptions about the event topology 
are tested. The decision is taken according to the obtained from the constrained fits 
with different topologies. 

Figure 1 shows the spectrum of the calculated energies for simulated and measured 
events passing all but the ^/s^ cut for 200 GeV e+e" centre-of-mass energy. The agree- 
ment between data and simulation is good for the high energies relevant to this analysis, 
while the peak around Mz appears to be slightly shifted in the simulation. A cut requir- 
ing the reconstructed centre-of-mass energy ^/s'^ to be greater than 0.9 • Ecm is applied 
to discard radiative return events (see Tabic 2). 

Two photon events are strongly suppressed by the cuts. Leptonic background was 
found to be negligible in this analysis. 
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neutral particle 
beiecLiuii 


E > 0.5 GeV 

ZU \ t/ ^ lUU 


charged particle 
selection 


0.4 GeV <p< 100 GeV 
Ap/p < 1.0 

measured track length > 30 cm 

H 1 Q.i' n r*(^ "I'o T T-' iTi T'/Ti TilfiTip ZL r'Tn 

LiloudliL/C \j\J ±.L . Ill / U/ UldliL/ \ t: L/lii 

distance to LP. in 2; < 10 cm 
20° < e < 160° 


Standard 

event 

selection 


A^ch>7 

30° < ^Thrust < 150° 

£;tot > o.50£;e„ 


ISR cuts 




WW cuts 


AT^h > 500 • S^in + 1.5 
A^ch < 42 


prompt photon selection 


- 10 GeV <Ew <E^ + h GeV 
11 GeV <E^< E^^/2 

= 25° 
E^ < 0.5 GeV 



Table 2: Selection of particles and events. E is the energy, p is the momentum, Ap 

its error, r the distance to the beam-axis, z the distance to the beam interaction point 
(LP.) along the beam-axis, and 6 the azimuthal and polar angles with respect to the 
beam, iVch the number of charged particles, 6'Thrust the polar angle of the Thrust axis 
with respect to the beam, E^ot the total energy carried by all particles, Ecm ther nominal 
Lep energy, y^i^ the reconstructed hadronic centre-of-mass energy, -Bmin is the minimal 
Jet Broadening (described in Section 2.3), Ej the energy of the detected photon, Ew the 
angular energy (see Equation 3), the opening angle of the photon isolation cone and 
Ea the maximum additional energy deposit within this cone. 
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Figure 1: Left: Reconstructed centre-of-mass energy for accepted data before \fs' cut 
compared to QCD and four-fermion simulation. Right: Simulation of four-fermion back- 
ground and QCD events in the iVch--Bmin plane. The upper two plots show the distribu- 
tions for semi-leptonic and fully hadronic WW events, respectively. The lines indicate 
the cut values chosen. 

Since the topological signatures of QCD four-jet events and four-fermion backgrounds 
such as hadronically decaying ZZ or WW events are similar, no highly efficient separation 
of these two classes of events is possible. Thus any four-fermion rejection implies a severe 
bias to the shape distributions of QCD events, which needs to be corrected with simula- 
tion. By applying a cut on an observable calculated from the narrow event hemisphere 
only (like -Bminj see Section 2.3), the bias to event shape observables mainly sensitive to 
the wide event hemisphere is reduced. The two dimensional cut in the A^ch^-Smin plane 
exploits the different correlation between these observables in QCD and four-fermion 
events (see Figure 1). Applying the two dimensional cut almost 90% of the four-fermion 
background can be suppressed. The remaining four-fermion contribution is estimated by 
Monte Carlo generators and subtracted from the measurement. The remaining detector 
and cut effects are unfolded with simulation. The influence of detector effects was studied 
by passing generated events (Jetset/Pythia [21] using the Delphi tuning described 
in [22]) through a full detector simulation (Delsim [16]). These simulated events are 
processed with the same reconstruction program and selection cuts as are the real data. 
In order to correct for cuts, detector, and ISR effects a bin-by-bin acceptance correction 
C, obtained from e"'"e~ — > Z/7 — > qq simulation, is applied to the data: 



gen,noISR 

Kfi) acc 



(2) 
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Figure 2: The energy distribution of the photon in the selected isolated photon events. 
The three classes (0,1,2) correspond to a mean centre-of-mass energy of 76, 66 and 45 GeV, 
respectively. The steps between the classes result from different tt^ rejection cuts. 

where /i(/i)gen,noisR represents the contents of bin i of the shape distribution / generated 
with the tuned generator. The calculation includes all stable particles. The subscript 
noISR indicates that only events without significant ISR {^/s — \fs' < 0.1 GeV) enter 
the distribution. /i(/i)acc represents the accepted distribution / as obtained with the full 
detector simulation. 

2.2 Data selection at hadronic centre-of-mass energies below 

In order to extend the available energy range below the Z-peak, events with reduced 
hadronic centre-of-mass energies due to hard photon radiation are selected from data 
taken at centre-of-mass energy of 91 GeV in 1992 through 95. The method requires 
an energetic isolated photon to be detected and is based on the hypothesis that such 
photons are emitted before or immediately after the Z/7 interaction [23] and do not 
interfere with the QCD fragmentation processes. The angular distribution of the initial 
state radiation is aligned along the direction of the beams, with the result, that most 
photons go undetected in the very forward region. In contrast, photons from final state 
radiation tend to be grouped along the direction of the final state partons and can be 
detected with better efficiency. As a result the selected events are dominated by final 
state radiation. 

Many of the photons convert into e"'"e~ pairs in the material before the calorimeter. 
These are reconstructed using the tracking detector information left by the e"^ and e~ 
particles. Only the reconstructed conversions before the TPC with electron and positron 
measurements in the TPC and conversions behind the TPC with electron and positron 
measurements in the OD or the HPC are used. 
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The largest part of the non-photonic background stems from tt s decaying into two 
photons. Due to the high granularity of the HPC the photon shower can be tested for 
a two photon hypothesis. This is done by two methods. The first tries to divide the 
cluster into two subclusters and reconstructs the invariant mass of the decayed particle. 
The second measures the asymmetry of the energy distribution in the ^^0-plane, as two 
photons generate a more elliptic cluster shape. The results of these two methods are 
combined in a single probabilistic variable. Since the angle between the two photons 
decreases with the energy, harder cuts had to be made for higher photon energies. 

In order to distinguish prompt photons from soft coUinear photons arising in the later 
stages of fragmentation and decays, hard cuts on the photon energy and the isolation 
with respect to other jets have to be applied. Isolation is defined by two criteria. In order 
to obtain photons at a large angle with respect to the final state particle a minimum angle 
of 20° to the jet axes is demanded. The jets are defined by the Durham cluster algorithm 
with Ucut = 0.06. The additional energy deposition within a 25° cone around the jet-axis 
had to be less than 0.5 GeV, which reduces background from 7r° decays. Electromagnetic 
punch-through entering the HCAL has been considered. 

To test the consistency of the measured photon energy, the following cross-check is 
performed: the event, cxluding the photon, is chistcrcd into two jets and the energy of the 
radiative photon is reconstructed from the angles between jets j,k and photon i through 
the following equation: 





sin 9jk 




1 sin 9ji 


+ 1 


sin 0ik\ 


+ 


sm6jk\ 



This reconstructed energy is required to coincide with the photon energy measured by 
the calorimeters in the range — 10 GeV < < + 5 GeV. 

The additional selection criteria for ISR and final state radiation (FSR) events are 
summarised in Table 2. The energy distribution of the final prompt photon candidates 
can be seen in Figure 2. From selected events the tagged photon is removed, and the 
event is boosted into the centre-of-mass frame of the hadronic system. The boost was 
determined by the measured photon. The events are summed up into three intervals in 
centre-of-mass energy. The mean value of each sample is taken as the nominal energy as 
calculated using the measured radiated photon and a correction is applied. 



2.3 Observables and their mass corrections 

The event shape observables used throughout this paper are calculated from the 
charged and neutral particles. 
The Thrust T [24] is defined by 

EN \ -> -> \ 

where pi is the momentum vector of particle i and ut is the Thrust axis, which maximizes 
the above expression. The observable Major is defined similarly, replacing nx with hm 
which is constrained to be perpendicular to nx- The C-parameter[25] is defined by the 
eigenvalues A of the infrared-safe linear momentum tensor Q^^ : 

I ^ i j 

Q^^' = v^iv i^i Elff C = 3(AiA2 + A2A3 + A3AO . (5) 
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Here pi denotes the i-component of Pk- Events can be divided into two hemispheres, 
positive and negative, by the plane perpendicular to the Thrust axis Ht- The so-called 
Jet Masses [26] are then given by 

M|=| J] Pfc) , M2 = max(M|,M!) , = Ml + 1\P . (6) 

y±Pfe-nT>0 J 

The Jet Broadenings [27] are defined by summing the transverse momenta of the particles: 
B± = ^ ^ \pk X utI , Einax = max(S_, B+) , Ssum ^ B_ + B+ . 

(7) 

The energy-energy correlation EEC [28] measures the correlation of the energy flow in an 
hadronic event. It is defined as 



VAcOSxf- ^ icos.-A-^ ^ 



dcosx N^^^^^^^^^ ^^jg ^cosx-=^ 



(8) 



Here Xij denotes the angle between the particles i and j. The jet cone energy fraction 
JCEF [29] integrates the energy within a conical shell of average half-angle x around the 
Thrust axis. It is defined as 



where Xt is the opening angle between a particle and the Thrust axis pointing in the 
direction of the light Jet Mass hemisphere: 

X. = cos . (10) 

Evidently the statistical correlation between the event shape observables when calcu- 
lated from the same data is very high. But since their QCD predictions are different (e.g. 
the relative size of the second order contributions) studying them provides an important 
cross-check for QCD. 

In the subsequent analysis QCD effects are calculated in the massless limit. This is 
an approximation which can lead to substantial deviations for some sensitive observables 
in certain cases and in the low energy limit in general. In order to reduce mass effects 
two approaches are applied: In [30] it is proposed that hadron masses are associated 
with corrections that are proportional to (In (5)^/(5 with A ~ 1.6, which can be of the 
same size as traditional power corrections. The mass induced power corrections can be 
separated into two classes, universal and non-universal, where the non-universal can be 
reduced by a redefinition of the observables. The Jet Masses in particular are subject to 
large mass corrections. To suppress the infiuence of the hadron masses, new observables 
were defined which for massless hadrons are identical to the standard M^/£'^;g-definition 
[30] . They are defined by replacing the four-momentum p in the standard formula: 

p={p,E) — > {pAp\) (p-scheme) (11) 
p = {p, E) — > {pE, E) (£;-scheme) (12) 
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Figure 3: Corrections due to mass effects: The full line shows the relative change of the 
distribution of the observable due to a change into the E-scheme as a function of the 
energy. The dashed hne shows the b-mass correction applied to the data. 
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with p being the unit vector in direction of p. 

The resulting observables have the same second order coefficients and the same power 
correction coefficient as standard definitions, because the theoretical calculations were 
performed for massless particles. Figure 3 shows, for shape observable means, the relative 
size of the corrections due to a change into the E-scheme as calculated by Pythia 6.1. 
M^/£'^jg and / E^^^ show the biggest changes (Other observables show in principle the 
same effect, but to a neghgible extent.). As a consequence the E definitions of M^/E^jg 
and M^/ E^^^ are used wherever possible. 

Besides the influence of final state hadron masses the influence of heavy b-hadron 
decays has to be taken into accoimt. In order to correct for the kinematic effects of b- 
hadron decays, samples of 100000 b-quark and 100000 light quark events were simulated 
using Pythia 6.1. for each energy used. For the calculation of the shape observables 
all stable particles were considered. A correction was then calculated as the ratio of the 
event shape distribution (or mean) for light quark to b events considering the energy 
dependence of the fraction of b events in e"'"e~ annihilation. The correction is applied to 
the data throughout the analysis. It is shown for several event shape means in Figure 3. 
While the overall effect on the observables is of the order of a few % at the Z-peak, it 
rises well above 5% at low energy and influences the evolution of the observables. 

2.4 Systematic uncertainties and definition of average values 

In order to estimate systematic uncertainties of the corrected data distributions and of 
the quantities derived from them, the event selection and the correction procedure were 
varied. For each variation the analysis was repeated. The individual deviations from the 
central result were added in quadrature and considered as the systematic uncertainty. 

The following variations were made in the event selection: the cut in the charged 
multiplicity was modifled by ±1 unit, the cut in the polar angle of the event Thrust 
axis was modified from 25° to 35°, and the cut on the observed visible energy was varied 
between 0.45 and 0.55. For the high energy data the y/s^ cut was lowered to 0.8. For data 
at centre-of-mass energies above the WW threshold the weight of Smin in the cut relation 
was lowered to 480 from 500 and the WW cross-section was increased conservatively by 
5%. 

When hadronisation corrections are included in the analysis the predictions of Ari- 
adne were used instead of the standard choice Pythia. Model parameters are as given 
in [22]. Additionally ±10% was taken of the hadronisation correction as the systematic 
uncertainty. For the kincmatically dominated b-hadron mass ±20% of the correction was 
taken as the systematic uncertainty. 

For all renormalisation scale dependences the scale / was varied between half and 
twice the central value. 

The fit results for the different observables are summarized by quoting two kinds of 
average values: the unweighted mean value with the R.M.S. as first error and the 
(error-)weighted mean value with the simple average of the individual statistical errors 
as the first error. The second error in both cases is the systematic uncertainty which has 
been propagated from the individual results. Quoting the R.M.S. indicates the size of 
theoretical uncertainties, while the average statistical error of the weighted mean value 
indicates the statistical significance of the results. Following from our definition, the sta- 
tistical error of the weighted mean value is bigger than some individual statistical errors. 
A treatment of statistical correlation has not been performed, since our errors are domi- 
nantly systematic. Furthermore the statistical correlation between e.g. the event shape 
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means is high (> 0.8), thus the gain in reducing the statistical error would be negligible. 
Using these highly correlated observables is useful for other reasons: It provides a cross- 
check for QCD and indicates theoretical uncertainties, since their perturbative expansions 
show different properties (e.g. a different size of the second order contribution). 

3 Comparison of experimental results to fragmenta- 
tion models 

Although among the oldest event shape measures, few data are available for the ob- 
servable Major. Therefore in the Figures 4 and 5 the Major distributions are shown for 
several energies between 45 GeV and 202 GeV compared to predictions of the Jetset and 
Ariadne Monte Carlo models. Except the lowest energy data at 45 GeV both simula- 
tions are almost indistinguishable. Data and simulation agree very well. The agreement 
between data and models is similarly good for other observables [31]. 

Figure 6 shows mean values in the energy range between 45 and 202 GeV for several 
observables, including the standard and the E definition of the Jet Masses. For compar- 
ison, results from Jetset simulations arc shown. Again good agreement between data 
and model is observed. The dotted line in Figure 6 represents the shape observable mean 
at the parton level ^. It is seen that the hadronisation correction, that is the difference be- 
tween hadron and parton level curves is smallest for the observables (Major), (M^/E^jg), 
and (-Bmax) ■ Also the slope of the parton and hadron level agrees best in these cases. On 
the other hand for the "subtracted" observables like (Oblateness), which have originally 
been constructed to compensate for hadronisation effects, show increased hadronisation 
corrections. In these cases the correction can have opposite sign to the other observables 
and sometimes even the sign of the slope of the energy evolution is opposite for parton 
and hadron level. 

The behaviour of the hadronisation correction indicates a clear preference for observ- 
ables such as (Major), (Ml/E^^^), and (-Bmax) which are mainly sensitive to the hard 
gluon radiation in the events. It should, however, be noted that in these cases some tech- 
nical problems may exist in the calculation of resummed theoretical predictions as soft 
gluon radiation may lead to a badly controlled exchange of the wide and narrow event 
hemispheres [30]. 

4 Power corrections to differential distributions 

When comparing event shape data to perturbative calculations in general corrections 
for the effects of non-perturbative hadronisation are applied. One approach to hadro- 
nisation corrections is the renormalon induced power correction model proposed in [9]. 
In this model the origin of non-perturbative effects is determined by Borel transforming 
the observables and fixing the singularities found on the real axis. For several differential 
distributions and in the simplest approximation this results in shifting the distribution 

^Parton and hadron level refer to the simulation of the hadronic final state. Parton level is before, hadron level after 
hadronisation has taken place. 
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Figure 4: Distributions of the observable Major for centre-of-mass energies of 45, 66, 76 
and 91.2 GeV compared to predictions of Jetset and Ariadne. In each plot the upper 
chart shows the size of the detector correction, defined as tt?^- The grey area indicates 
the distribution of non-radiative background events which has been subtracted. 
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Figure 5: Distributions of the observable Major for centre-of-mass energies of 189, 192, 
200 and 202 GeV compared to predictions of Jetset and Ariadne. In each plot the 
upper chart shows the size of the detector correction, defined as The grey area 

indicates the distribution of WW and ZZ background events which has been subtracted. 
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Figure 6: Event shape means for different observables in comparison to Pythia 6.1 
predictions. Tfie full line shows the hadron level, the dashed hne the parton level. 
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T>f[f) of the observable /: 

Vf{f)^VY\f-V-Cf) , 



In 



Er 



K 



+ 7^ + 1 

1^1 Po 



(13) 



K = (67/18 - 7rV6)CA - 5nf/9 



Cp and Ca are the QCD colour factors, ckq is a non-perturbative parameter accounting 
for the contributions to the event shape below an infrared matching scale /ij. The Milan 
factor M. for three active flavours in the non-perturbative region is 1.49 [32]. 

Since the derivation of the coefficient c/ is based on resummation, there are only pre- 
dictions for exponentiating observables available, c/ is an observable-dependent constant, 
which is identical to the coefficient in the predictions for event shape observable means. 
These are: 



Observable 


1 -T 


C 




S / VIS 




2 


37r 


1 


2 



In order to show all formulae in a coherent fashion, we use the definition of [1] for the 
coefficients of the /3-function: 



/3o = 



33 - 2nf 



153 - 19nf 



132 = 2857 



5033 



-rii + 



325 



-rir 



3 ' ' 3 ' ' 9 27 

4.1 Power corrections for the Jet Broadenings ^max and B. 



(14) 



sum 



Unlike the former observables, the Jet Broadenings cannot be sufficiently described by 
simple shifts, as the shift becomes a function of the Jet Broadening. Early predictions 
neglected the recoil of the quark due to the gluon emission, which proved to be an 
important effect. Improved calculations [33] take this mismatch into account. For the 
wide Jet Broadening -Bmax the correction coefficient has the form 



CB^ax(^ma.) = 2 (Vlni^„.ax + ^70 - 2 - pi^R!) + x(^') + V(l + ^') + V'(l)) , (15) 
where 



= 2C^ ^-^^— ^^ (In^-^ - ^) , = -0.6137 , 

TT 4 



X(a) = - ([A(a)r - 1) , 



For the total Jet Broadening -Bgum the correction factor is 



CB_(Ssum) = 2 (2cB.3.(^sum) + 2[V'(1 + 27^') - V^(l + U')] + H{B-')) (16) 



with 



H{x) = £^eMn^)-nz)}Y^ 



^(l + 27^') 25sum 



+ 2n' + n'{z)) 



e^^A(7^') ' 
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with 7e = 0.5772 being the Euler constant and Zq given by the position of the Landau 
Pole of the two loop radiator TZ: 



n{x) 



4Cp 



c- 



In 1 - 



In a; 



+ lnx 



27r 



PoOis (l 



27r 



) 



4.2 Power corrections for the Energy Energy Correlation EEC 

The power corrections for the EEC have been calculated in [34] . Unlike for the other 
observables, there is no simple factorisation of the perturbative and non-perturbative 

components possible. Instead the non-pcrturbative coefficient is a part of the radiator 
function. The dominating non-pcrturbative part is based on the quark-gluon radiation 
and in the limit of large angles has a Q~'^ like behaviour, with 7 ~ 0.32 for Uf = 5. As a 
result one gets 



;i + tan2|)3Q 



d cos X 4 2 

where the perturbative radiator is defined by 



/ Jo(6gtan§)e-^'''^'(''^-5%£(l-26A) , (17) 

Jo 2 



etc 



-(ln(l -£)+£) 



3^ 
' Stt 



ln(l -£) + 



Pi 



2nPo V2 



in(i -e) e 



1-i 



+ 



1-e 



Po — in 

'^"27r 2 



The linear non-pcrturbative correction —2bX stems from the correlation between quarks 
and soft gluons, where A characterises the non-perturbative interaction at small momen- 
tum scales: 



A 



4Cf 



(18) 



The coefficient C{as) has the form 



11 tt' 



27r 



The radiator has its own non-perturbative component a: 



Cf 2 



PT{ 



(19) 



The treatment of this non-perturbative component is unclear, as it is much weaker than 
the linear part and it is missing the quadratic part from the quark-gluon correlation. The 
^pqil^i) normalisation factors: 



a, 



PT 



a. 



Po_ 
27r 



Q 



, ^ K q+1 
In — -I h 

fjLi Po P + 1 



a. 



(20) 



As a result one has to deal with three non-perturbative parameters: ao,ai and ai^i, 
where ao is equivalent to the non-pcrturbative parameter of the other observables. A 
phenomenological prediction of ai is given in [34] from DIS experiments and the theo- 
retical predictions are ai — 0.45 and ai^i — 0.55. Since the formula for the EEC is a 
large angle approximation it is important to choose a proper fit region. Below 120° the 
prediction becomes invalid, while above 170° the infiuence of higher order logarithmic 
terms can no longer be neglected. 
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4.3 Comparison with data 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.05 0.1 0.15 0.2 0.25 03 



Figure 7: Jet Broadening distributions as measured by DELPHI for centre-of-mass en- 
ergies between 45 and 202 GeV. The full line indicates the power model fit in the fit 
range, while the dotted line shows the extrapolation beyond the fit range. The dashed 
line shows the result after subtraction of the power correction. 

The derivation of the power correction predictions in general relies on the resummation 
of logarithmically divergent terms. The validity of these predictions is thus limited to 
a kinematical region close to the two jet regime. This has been taken into account 
when choosing the fit intervals indicated in Table 3. Additionally it was required that 
the corrections applied to the data as well as the fit results obtained were stable. The 
experimental systematics have been determined as discussed in Section 2.4. In addition, 
changes of the fit ranges were apphed as tabulated in Table 3. All systematic studies enter 
into the specified systematic uncertainty. Also the so called R matching was applied for 
the perturbative prediction instead of the standard logR matching [35]. It is notable that 
the change of the matching scheme has a significant influence on the size of the power 
corrections. 

Examples of fits compared to the data on the Jet Broadenings, 1 — Thrust and the 

EEC, are shown in Figures 7 and 8. Results of the fitted parameters for all observables 
determined from DELPHI data are given in Table 4. The expected correlation of the fit 
parameters ag and cto is displayed in Figure 9. The /ndf of the fit is acceptable when 
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Figure 8: EEC and 1 — Thrust distributions as measured by DELPHI for centre-of-mass 
energies between 45 and 202 GeV. The full line indicates the power model fit in fit range, 
while the dotted line shows the extrapolation beyond the fit range. The dashed line shows 
the result after subtraction of the power correction. 

systematic uncertainties are included. The as values tend to be rather low compared 
to the world average value of as{Mz) = 0.1181 ± 0.002 [1] for most observables. The 
non-perturbative parameter is higher for the classical M^/E'^^^ as expected due to the 
influence of hadron mass effects. For the other observables the results for ao agree within 
a relative uncertainty of about 20%. 

The result of the fit of the EEC is shown in Figure 8 (left). The influence of the non- 
perturbative parameters ai and ai^i is found to be much smaller than the influence of ao- 
The precision of the experimental data is insufficient to determine ai^i. A three parameter 
fit neglecting yielded: a, = 0.1173 ±0.0021 ±0.0008, ao = 0.478 ± 0.046 ± 0.017 and 
ai = 0.005 ± 0.026 ± 0.025 with a x^/ndf = 52.7/90. Since this fit indicates that the 
non-perturbative part of the radiator ai can be neglected, an additional two parameter 
flt was performed, resulting in: as = 0.1171 ±0.0018±0.0004, ao = 0.483±0.040±0.011, 
with /ndf = 53/91. The ao value for the EEC is consistent with the values determined 
from the other observables (see Table 4) . 
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Observable 


lower 


upper 


1 — 1 


n no ( r\ ^o^ 
U.Uo (^U.Uzj 


0.2 


(0.24) 


-'-'max 


U.Uo (U.U4j 


U.i4 


(U.lbj 


D 

-'-'sum 


U.Uy (U.Uoj 


n 1 T 
U.i / 


/ c\ 1 ^^ 

(u.iyj 


/\/f2 / I7i2 


U.Uo (^U.Uzj 


U.iz 


(U.14j 


/\/f2 / 7712 


U.Uo (^U.Uzj 


0.12 


(0.14) 




0.03 (0.02) 


0.12 


(0.14) 




0.03 (0.02) 


0.12 


(0.14) 


M!^/El, 


0.03 (0.02) 


0.12 


(0.14) 


Ml/El, 


0.03 (0.02) 


0.12 


(0.14) 


C-parameter 


0.2 (0.16) 


0.68 


(0.72) 



Table 3: Fit intervals used for the fit of power corrections to event shape distributions. 
The variations of the fit interval used for systematic studies are shown in brackets. 

4.4 Measurement of the non— perturbative shift from the Su- 
dakov Shoulder 

So far all predictions apply near to the two jet region. For some observables, however, 
the predictions seem to hold even in the far three jet region. In order to provide evidence 
for this observation the three jet limit of several observables was studied. 

Despite the fact that QCD event shape observables are constructed so as to be infrared- 
and coUinear safe, there can still be infinities at accessible points in phase space. The 
finiteness is only restored after the resummation of divergent terms to all orders. The 
resulting structure is called a Sudakov shoulder [36]. The most common case is the 
phase space boundary of the three jet region, which introduces a visible edge into the 
distributions. The position can be calculated and is, for example, 2/3 for 1-Thrust and 
3/4 for the C-parameter. A simple inspection of the distributions and the corresponding 
model curves shows that the shoulder is typically shifted to higher values, and that power 
corrections describe the shift rather well. 

The shift for the C parameter can be measured by fitting the slope of the logarithmic 
distribution on both sides of the shoulder. The intersection of these fits is a good ap- 
proximation to the shoulder position. The result of the fit can be seen in Figure 10, the 
fitted position of the shoulder is at C= 0.794 ± Q.Ql&stat i O.OOls^s corresponding to a 
shift of +0.044 ± 0.016stat ± 0.001^^^, with respect to the nominal position. Using a value 
of cc = 37r for the C parameter, a value of = 0.476 ± 0.097stat =t O.OOlbgys is obtained 
from Equation 13. This result is well consistent with the result obtained from the fit of 
the overall distribution ao = 0.502 ± 0.005stat ± 0.047^^5 suggesting a constant shift over 
the whole three jet region in the case of the C parameter. 
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Observable 




ao{iJ,i ^ 2GeV) 


X /ndf 


1-T 


0.1154±0.0002±0.0017 


+0.0004 


0.543+0.002+0.014 +0.013 


291/180 


Eiaax 


0.1009±0.0003±0.0016 


-0.0018 


0.571+0.005+0.031 +0.021 


106/90 


-^sum 


0.1139±0.0006±0.0015 


-0.0035 


0.465+0.005+0.013 +0.008 


88/75 


Tl t'? / 7-l9 


0.1076±0.0001±0.0013 


+0.0003 


0.872+0.000+0.026 +0.005 
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^E/^vis 


0.1056±0.0003±0.0006 


+0.0001 


0.692+0.007+0.010 +0.010 


120/90 


71 /f 2 / I7'2 

^^hp/-^vis 


0. 1055±0.0004±0.0010 


+0.0001 


0.615+0.009+0.022 +0.010 


130/90 


M hi ■ 

SI VIS 


0.1190±0.0004±0.0030 


+0.0001 


0.734+0.004+0.034 +0.009 


66/45 


/\/r2 / 1712 


0.1166±0.0004±0.0028 


+0.0002 


n t^cQ_i_n nn/i_i_n nov i n nnv 
U.boo±U.UU4±U.Uz / +U.UU/ 


OU/4o 


Ml/El, 


0.1156±0.0005±0.0010 


+0.0001 


0.536+0.005+0.010 +0.008 


54/45 


C-Paramcter 


0.1097±0.0004±0.0032 


-0.0008 


0.502+0.005+0.047 +0.021 


191/180 


weighted mean 


0.1078±0.0005±0.0013 


-0.0012 


0.546+0.005+0.022 +0.013 




unweighted mean 


0.1110±0.0055±0.0007 


-0.0008 


0.559+0.073+0.009 +0.013 




EEC 


0.1171±0.0018±0.0004 


0.483+0.040+0.011 


53/90 



Table 4: Determination of ag and cto from a fit to event shape distributions. Only Delphi 
measurements are included in the fit. The first error is the statistical uncertainty from 
the fit, the second one is the systematic uncertainty, the third the difference with respect 
to the R matching scheme. Only E-definition Jet Masses have been taken for the means. 
For the definition of the mean values see section 2.4. 



5 Power corrections in event shape means 

The mean values of event shape variables are defined as: 

(7tot J da A^evt. 

We have calculated them from the detector corrected and binned distributions. Hence 
they are fully inclusive quantities depending on a single energy scale only, and are well 

suited for low statistics analyses as the statistical uncertainty is minimised by using all 
events. Though the characteristics of the event shape obscrvablcs may differ in specific 
regions of the value of the observable, global properties can be assessed from the energy 
dependence of the mean value. 



5.1 The Dokshitzer and Webber ansatz 

The analytical power ansatz [12,37] including the Milan factor [10,11] is used to de- 
termine as from mean event shapes. This ansatz provides an additive non-perturbative 
term (/pow) to the perturbative 0{a^) QCD prediction (/pert), 

(/) = (/pert) + (/pow) , (21) 

where the 2nd order perturbative prediction can be written as 

y„„.H,Mi),(,.,„,„£.,,)(^^)\ 



21 



0.8 



0.6 



0.4 



0.2 





X 




DELPHI 




□ 




•1-T 








, C-Parameter 

















0.09 0.1 0.11 0.12 0.13 0.14 



Figure 9: Results of the fits to shape distributions in the a^-ao plane. The band indicates 
the world average of a^. 



A and B are known coefficients [38-40] and fi is the renormalisation scale. The power 
correction is given by 

(/pow) = c/P , (23) 

where V is as defined in Equation 13. 

The observable-dependent coefficient is identical for shapes and means. In the case 
of the jet broadenings Cf cannot be described as a constant. Here the non-perturbative 
contribution is proportional to l/{Q^/ces{Q)) [33]: 

I T^yCB , 3 ISqCb \ , . 

= cb — , +-- -rpr- + ^0 , (24) 

where cb is 1/2 in the case of (-Bmax) and 1 for (-Bgum), ^70 = —0.6137056. 

In the following analysis the infrared matching scale /ij was set to 2 GeV, as suggested 
in [12], and the renormalisation scale /i was set to -Ecm- 

A combined fit of as and ao to a large set of measurements^ at different energies [41] 
has been performed. In the calculation, statistical and systematic uncertainties were 
considered. For E^m > Mz, only Delphi measurements were included in the fit. Fig- 
ure 11 (left) shows the measured mean values of (1 — T), {M^^/ E^^^) (standard- and E 

■^We concentrated on up-to-date results. 



22 




Figure 10: Determination of the shift of the Sudakov shoulder in comparison to predic- 
tions. The vertical positions are arbitrarily scaled. The fit range (straight line) and the 
extrapolation to the intersection point (dotted line) are shown. The vertical line denotes 
the three jet limit of 0.75 . 

definition), (C), (-Bmax) and (-Bsum) as a function of the centre-of-mass energy together 
with the results of the fit. The fit values of and are summarised in Table 5 and 
displayed in the a^-ao plane in Figure 11 (right). The systematic uncertainty was ob- 
tained as described in Section 2.4. In addition /xi was varied from 1 GeV to 3 GeV. Both 
uncertainties were added in quadrature. 

The Us values obtained from these fits are consistent with each other and in good 
agreement with the world average as{Mi) = 0.118 ± 0.002 [1]. The extracted ao values 
are around 0.5 as expected in [11,37]. However, the predicted universality (e.g. the 
observable independence) is satisfied on a 25% level only. This problem remains even 
though only the Jet Mass in the E-definition, {Ml^/E^^^), is considered which avoids 
a strong additional energy dependence due to the influence of hadron masses. The as 
values are higher and the ao values are lower than the corresponding results from event 
shape distributions (compare Figures 9 and 11). 
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Observable 


0(o[iJ>i — ^ (jew ) 


«s(Mz) 


xVndf 


(l-T) 


0.491±0.01d±0.009 





1241±0.0015±0.0031 


26 


f / A 1 

5/41 


(C- Parameter) 


r\ A A A _i_r\ r\or\_i_r\ r\r\o 





1222±0.0020±0.0030 


11 


6/23 


/ A/r2 / T7<2 \ 
(^h/^vis) 


0.601±0.058±0.012 





1177±0.0030±0.0018 


14 


1/27 




0.300±0.222±0.127 





1185±0.0104±0.0057 


10 


1/15 


/ ]\j2 / rn2 \ 


u.ooy±(j.zzy±(j.izy 





1197±0.0107±0.0058 


9 


5/15 


/ ]\/r2 / TP2 \ 


U.544±U.ioU±U.Uyo 





1335±0.0118±0.0074 


1 


// io 


/ A/f 2 / T?2 \ 


U.o ( oitU. iooitU.Uo4 





1288±0.0104±0.0067 


8 


7/15 




0.409±0.143±0.086 





1304±0.0107±0.0069 


8 


2/15 


(-^max) 


0.438±0.041±0.027 





1167±0.0018±0.0007 


10 


1/23 


{Bsnm) 


0.463±0.032±0.009 





1174±0.0021±0.0020 


8 


8/23 


weighted mean 


0.468±0.080±0.008 


0.1207±0.0048±0.0026 






unweighted mean 


0.431±0.048±0.039 


0.1217±0.0046±0.0030 







Table 5: Determination of cto and ag from a fit to a large set of event shape mean values 
measured from different experiments [41]. For Ecm > Mz only Delphi measurements are 
included in the fit. The first error is the statistical uncertainty from the fit, the second 
one is the systematic uncertainty. For the mean values only the E-definition Jet Masses 
have been used. For the definition of the mean values see section 2.4. 

5.2 Simple power corrections 

Power corrections in the Dokshitzer- Webber framework can only be calculated for the 
set of exponentiating observables. The experimental evidence for corrections which show 
1/Q behaviour is however not restricted to this type of observables. The tube model 
indicates the existence of power corrections on simple phase space assumptions. 

In order to determine approximate power corrections for all observables measured the 
"simple power correction" ansatz is used. Here an additional power term (fpow) = Ci/Q is 
added to the 0{al) perturbative expansion of the observable, with Ci being an observable 
dependent, unknown constant. The disadvantage of this simple ansatz is that a double 
counting in the infrared region of the observables is not corrected for as in the Dokshitzer- 
Webber approach. 

Two different types of fits were performed to the DELPHI data using this simple model. 
Firstly, in order to investigate whether this simple model yields sensible values for ag at 
all, both parameters, Aj^g and Ci, were left free in the fit. The results obtained are given 
in Table 6. The average value of these as results and the corresponding R.M.S. obtained 
only from the fully inclusive observables are ag— 0.1250± 0.0058 for the unweighted mean 
and as— 0.1250 ± 0.0054 for the weighted mean. The reasonable values for ag as well as 
the acceptable /ndf of 6/7 of the averaging support the approximate validity of this 
simple power correction model. Secondly, in order to get comparable estimates of the 
size of the power correction for the different observables, A-ms = 0.250 GeV was chosen, 
leaving only Ci as a free parameter. The fitted values of Ci are contained in Table 7, 
where the total experimental uncertainty for Ci is given. 

The ratio of the power model parameter normalised to the first order perturbative 
coefficient, Ci/A, is plotted against the ratio of second to first order perturbative co- 
efficient B/A in Figure 14. The normalisation to A is made in order to make the ob- 
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Observable 


Ci 


as{Mz) 


xVndf 


(l-T) 
\^iviajor/ 


0.483±0.085±0.048 
u.4uo±u.y / y±u.oyo 


0.1312±0.0019±0.0032 
nil fifi-un nn7i -un nn/i q 

U.iiDDztU.UU/ i±U.UU4o 


26.3/41 

S /I 

y.o/ io 


(C-Paraiiietor) 


2.212±0.169±0.219 


0.1278±0.0021±0.0032 


11.0/23 


/ ;i/f2 / zr'2 \ 

(iWhE/^vis) 
V'^^sEZ-^vis/ 


U.oUzdiU.iiydiU.Uz / 

0.736±0.689±0.408 

U.ioU±U.4o4±U.z4D 

n Q/11 -i-n fii 7-i-n '?7'? 

U.O^ilCU.Oi ( ICU.O ( 


n 1 01 n_Ln nnQQ_Ln nni q 

0.1375±0.0132±0.0081 

U.iziorhU.UiiozhU.UUDi 

n 1 Q/io-un ni 1 Q-i-n nn7c; 
u. io^izicu.ui iyicu.uu ( o 


io.i/z/ 
7.1/15 
y.o/ io 

0. i/ iO 


( -Smax ) 

/ R \ 
\J-^sum/ 


0.241±0.075±0.018 
u.oyoinu. ioymu.uou 


0.1203±0.0016±0.0009 
n lOQK-i-n nnis-i-n nnoo 

U.IZOUIlU.UUIOIlU.UUZZ 


9.2/23 
( .0/ zo 


(EEC7oo-iio°) 

(EEC3oo-l50°) 


0.285±0.637±0.541 
0.022±0.470±0.691 


0.1307±0.0102±0.0088 
0.1395±0.0051±0.0082 


18.8/15 
41.2/15 


(JCEFiioo_i6o°) 


0.011±0.676±0.954 


0.1201±0.0049±0.0070 


31.4/15 


weighted mean 
unweighted mean 




0.1250±0.0054±0.0024 
0.1250±0.0058±0.0032 





Table 6: Determination of Ci and ag from a fit to a large set of measurements of different 
experiments [41]. For > Mz only Delphi measurements are included in the fit. 

The first quoted error is the uncertainty from the fit, the second one is the systematic 
uncertainty. In calculating the mean values, the Jet masses using standard definitions, 
EEC and JCEF have been omitted. For the definition of the mean values see section 2.4. 



Observable 


Ci 


(l-T) 
(Major) 


1.093±0.016±0.289 
0.390±0.020±0.170 


(C-Paraiiieter) 


1.330±0.0(33±0.713 


(Ml/El,^) 

(M?e/^v\) 


0.649±0.023±0.076 
1.839±0.051±0.136 
0.327±0.018it0.072 
1.265±0.033±0.136 


(-^max) 
{Esum) 


0.406±0.008±0.046 
1.199±0.012±0.190 


(EEC7oo-llo°) 
(EEC3o°_i50°) 


1.177±0.021±0.163 
2.066±0.015±0.227 


(JCEFiioo-i6o°) 


0.484±0.026±0.223 



Table 7: Determination of Ci with a fixed = 0.250 GeV from a fit to a set of 

measurements of different experiments [41]. For E^m > Mz only Delphi measurements 
are included in the fit. The first uncertainty is the uncertainty from the fit, the second 
one is the systematic uncertainty. For the definition of the mean values see section 2.4. 
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Figure 11: Left: Measured mean values of (1 — T), (M^/i?^jg), (C), (-Bsum) and (-Bmax) 
as a function of the centre-of-mass energy. For clarity some of the high energy data have 
been merged. The solid lines present the results of the fits with Equations (21-23), the 
dotted lines show the perturbative part only. Right: Results of the Dokshitzer- Webber 
fits in the Os-ao plane. The band indicates the world average of a^. 

servables directly comparable. Experimentally a clear correlation between the genuine 
non-perturbative parameter Ci and the purely perturbative parameter B is observed. 
This strong correlation indicates that the term Ci/Q should not be interpreted as purely 
non-perturbative. 
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Interpretation of event shape means using RGI 
perturbation theory 



Today the MS scheme is commonly used for the representation of perturbative cal- 
culations of physical observables. In consequence predictions for power corrections have 
also been given in this scheme. Nonetheless the MS scheme is only one of an infinite set 
of equally well suited schemes. 

A previous Delphi analysis using experimentally optimised scales for the determina- 
tion of as from event shape distributions [42] has shown that the experimentally optimised 
scales of different observables are correlated with the corresponding effective charge, ECH 
[43], or principle of minimal sensitivity, PMS [44], scales. The renormalisation group in- 
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variant (RGI) approach [5,7] uses the same central equations as the method of effective 
charge [43], however the motivation and philosophy differ. The derivation of the RGI 
method makes no reference to any renormalisation scheme whatsoever. In the follow- 
ing sections the RGI predictions [5,45,46] for fully inclusive shape observable means are 
explored. 



6.1 Theoretical background of RGI 

Instead of expanding an observable R into a perturbative series in as{Q), the starting 
point of the RGI method is a /3-function like expansion of ^ in R: 

Q— = -bR^l + p,R + P2R^ + ...) = bp{R) . (25) 
dQ 

Note that here R — 2{f)/A is normalised such that the perturbative expansion in ag/Tr 
begins with unit coefficient. It can be shown [5,43], that the coefficients pi are scheme 
invariant and that the scale dependence cancels out completely. The pi can be calculated 
from the coefficients ri,r2 of the perturbative expansion: 

Pi = ci , P2 = C2 + r2 - rici - rl , (26) 

,_(3o _ /3i _ /32 _ B _ C 

The coefficients A and B are defined in Equation 22, C is the corresponding third order 
coefficient. 

The solution of Equation 25 is equivalent to the well known implicit equation for a^: 

bln-^ = --p^\n(l + ^]+ rdx(^ + ^-^ 7^ ■ (27) 

Ar R \ piRj Jo \pix) x^il + pix)J ^ ' 

Here A.r is an R-specific scale parameter. In next-to-leading (NLO) order as the integrand 
vanishes and the solution of this equation is identical to a scale which sets the NLO 
contribution of the perturbative series to zero. Using the so called Celmaster Gonzalves 
equation [47] , can be converted into Ajjg without any loss of precision: 



A 



e-/^ -i . (28) 



In a study of event shape observables a test of the validity of RGI perturbation theory 
is currently limited to NLO. It should therefore be verified that higher order corrections 
are small and that the NLO /3-function (e.g. the inclusion up to the pi term in Equation 
25) is a sufficient approximation. A check on higher order contributions is implied by a 
consistency check of the values measured from different observables. 

It is important to note that the above derivation only holds for observables that depend 
on one single energy scale, such as fully inclusive observables. Selections or cuts in the 
observable introduce additional scales that have to be included into equation 25. Thus 
it is not to be expected that the aforementioned simple form of RGI is valid, say, for 
Ucut dependent jet rates or bins of event shape distributions. The application of the RGI 
method to ranges of the EEC or JCEF as performed in the following sections is thus not 
fully justified, except by the success of the comparison to data. The RGI method may 
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also to some extent apply here as the intervals chosen are rather wide and represent an 
important fraction of the events. Moreover it should be noted that the total integrals 
over the EEC or JCEF are normalised to 2 or 1, respectively. RGI can also not be 
unambiguously calculated for observables like Mj — — M^, as Mi — min(M^, M^) is 
only known to leading order. 

It is possible to include power corrections into RGI. In [45] it is shown that the exis- 
tence of non-perturbative corrections leads to a predictable asymptotic behaviour of the 
renormalisation group equation. This can be included into the equation as a modification 
of the p function [48] : 

p(x) ^ p(x) - ^x-^^/'e-'/'^ , (29) 

where Kq is a free, unknown parameter that determines the size of the non-perturbative 
correction. The correction is approximately equal to a simple Ci/Q power correction in 
the MS scheme with [48]: 

= -i^o • e-/^ • j -Ams ■ (30) 

As RGI pertubation theory and the ECH method are based on the same basic equa- 
tion, the choice of the ECH renormalisation scheme is implicit in RGI perturbation theory. 
Therefore it may be controversial whether measurements of performed using RGI are 
renormalisation scale or scheme independent. However, with respect to the (3 function 
(e.g. its leading coefficients, /5o and [3i) the situation is different. Their measurement 
based on Equation 25 is free of any scheme ambiguity since this relation holds in any 
renormalisation scheme. Moreover /3o and /3i are renormalisation scheme invariant quan- 
tities. Dhar and Gupta summarize their discussion with the following words: "(...) we 
have shown that in a renormalizable massless field theory with a single dimensionless 
coupling constant, only the derivative p{R) of a physical quantity R with respect to an 
external scale is well defined and unambiguously calculable" [5] . 

6.2 Comparing RGI with power corrections to data 

Using the same data as in the case of simple power corrections, a combined fit of Ar 
and Kq is performed to the RGI with power correction theory using Equation 29. A 
correction is applied for the influence of the b-mass but not for further hadronisation 
effects. RGI plus power correction describes the behaviour of the data well for all observ- 
ables considered including EEC and JCEF and leads to a consistent result for ag- The 
fit results are given in Table 8 and Figure 12. The uncertainties for the jet masses in 
E-definition, Major as well as for JCEF and EEC are large as no data from low energy 
experiments are available in these cases. The result for (1 — T) agrees reasonably with the 
comparable analysis presented in [48]. Note that in [48] no b-mass correction was applied. 
Small differences are understood as being due to the b-mass correction and the differing 
choice of input data. Combining the results for the fully inclusive observables yields 
as= 0.1179 ±0.0020 ±0.0013 for the unweighted mean and a,= 0.1184 ± 0.0031 ± 0.0035 
for the weighted mean with a /ndf of 7.0/5. The spread of the otg results from the fit 
of RGI with power corrections is only half of the size as for the simple power correction 
case. 

The fact that the results for Kq are small (most are compatible with zero) is surprising. 
The smallness of the Kq results in combination with the consistent values from all 9 
different observables casts doubt on the interpretation of the measured (fpow) contribution 
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Figure 12: Results for a combined fit of and the power correction parameter Kq using 
Equation 29. Left: The results for the non-perturbative parameter Kq. Right: The 
results for as{Mz) deduced from Ar. The straight line shows the unweighted mean, the 
shaded band the R.M.S. obtained from the fully inclusive observables only. 



(see Equation 21) as a genuine non-perturbative term. In view of the observed consistency 
one is led to presume a better approximation of the data by RGI perturbation theory and 
the MS scheme appears as an unfortunate choice. This conjecture implies that the power 
terms measured for shape observable means in the previous section mainly parameterise 
terms which can in principle be calculated perturbatively. 

To elucidate this conjecture further, Thrust is taken as example. Here for a fixed 
as=0.118, Ci/Q = (1.108 ± 0.017) GeV/Q is found as the simple power correction in 
the MS scheme. Using Equation 30 this translates into Kq = —0.32 ± 0.006. In the 
RGI plus power term approach Kq = —0.005 ± 0.008 is observed, however. Contrary 
to the simple power model or the ansatz of Equation 23, here the power contribution 
is insignificant. From the influence of hadron masses, a small kinematical dependence 
similar to an inverse power law is to be expected. However, the absolute size of this mass 
correction cannot yet be safely calculated [30]. 



6.3 Comparing pure RGI to data 

In light of the negligible power term (cx Kq) observed in the previous section, a com- 
parison of the data and the pure RGI perturbation theory appears interesting. No hadro- 
nisation correction is included and Ar is the only free parameter of the theory. This 
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Observable 


Ko 




X^ndf 


(1 -T) 

(Major) 


— U.UU5±U.(J(Jo±(J.(JUo 
-0.126±0.302±0.197 


r» 1 1 n /I _i_n nnnn_i_n nnno 

U.1194±U.UUU9±(J.(J(J(Jo 
0.1143±0.0052±0.0033 


31.U/39 
9.8/13 


(C-Parameter) 


— (j.(j4(j±(j.(jiD±(j.(j(jy 


U.1175±(J.(J(J12±(J.UU1(J 


ll.U/21 


/ T\/r2 / 7712 \ 

{^hE/Eyis) 

{M^e/ED 


-0.112±0.028±0.006 
— U.(Jo4±(J.(J7z±(J.(J35 

— U.U /4±U.zoo±U.i / z 
-0.018±0.054±0.061 


0.1219±0.0014±0.0005 
U.1249±U.(J(Jd5±(J.(J(J35 

U.iiOo±U.UUoo±U.UUoU 
0.1205±0.0040±0.0029 


19.6/25 
7.I/I0 

n Q /I Q 

y.o/ io 
11.0/13 


1 R \ 
\-'-'max/ 

{Bsum} 


— U.ouz±u. iUoitU.Uzi 
-0.055±0.022±0.004 


nil ofi-i-n nni n-un nnnc; 
u.iiyo±u.uuiu±u.uuuo 

0.1169±0.0009±0.0005 


0.0/zi 
7.7/21 


(EEC70°-110°) 
(EEC3oo-150°) 


0.008±0.095±0.075 
0.045±0.036±0.054 


0.1174±0.0065±0.0051 
0.1236±0.0033±0.0047 


18.8/13 
40.8/13 


(JCEFiio°-i6o°) 


0.024±0.171±0.254 


0.1160±0.0043±0.0064 


31.4/13 


weighted mean 
unweighted mean 


-0.018±0.114±0.068 
-0.097±0.114±0.057 


0.1184±0.0031±0.0035 
0.1179±0.0020±0.0013 





Table 8: Results for a RGI plus non-pcrturbative parameter fit to a large set of mea- 
surements of different experiments [41]. For > Mi only DELPHI measurements are 
included in the fit. The first uncertainty is the statistical uncertainty from the fit, the 
second one is the systematic uncertainty. For the mean values only the E-definition Jet 
Masses have been used and both EEC and JCEF have been omitted. For the definition 
of the mean values see section 2.4. 

implies that can be precisely determined from each individual shape observable mean 
measurement allowing for stringent tests of RGI by comparing the energy dependence 
of individual observables or by comparing the A^g^ or as{Mz) values obtained from dif- 
ferent observables. Note that these parameters can be deduced from Ar without loss of 
precision (see Equation 28). 

Observed deviations between different energy points of an observable should mainly 
be due to hadronisation or mass corrections. From the above section these corrections 
arc ah^cady known to be small (~ few %). Differences in as values between different 
observables represent different power corrections as well as missing higher order terms of 
RGI perturbation theory. Clearly the ag value for each observable will be dominated by 
the precise measurements at the Z pole. 

Figure 13 compares the energy dependence of the fully inclusive event shape means 
considered in this analysis to the RGI prediction. The data shown are the same as ex- 
plained in the simple power correction section. The energy dependence of all observables 
is well represented by RGI theory. The 'X^ /ndf of the fits is acceptable for all observables 
(see Table 9) and similar to the RGI plus power correction case. The observed values of 
as of the different observables agree reasonably well, with the exception of the standard 
definition of M^/E^^^ or / E^^^ (not shown here). For these observables mass terms are 
known to be important and reflect into higher as values compared to the results using 
the E-scheme. 

The RGI as results for the different observables are also shown in Figure 15 compared 
to results of RGI with power corrections and values obtained in MS with simple 
power corrections as well as with Monte Carlo hadronisation corrections. The last have 
been calculated as usual by comparing the hadron and parton level result in a Monte 
Carlo fragmentation model (P ythia) . It is evident that RGI leads to a far better con- 
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Figure 13: Comparison of the data on event shape means with the prediction of pure 
RGI perturbation theory (full line). For clarity some of the high energy data have been 
merged. The dashed line represents the MS expectation with the same ^^(M^). 

sistency between different observables than the simple power corrections or the Monte 
Carlo hadronisation corrections. In comparison to the RGI plus power correction ansatz 
the spread of the values is not significantly increased. The unweighted mean and 
R.M.S. spread of the values from the fully inclusive means shown in Figure 15 is 
as= 0.1201 ± 0.0020, using the weighted mean yields 0.1195 ± 0.0002. The average 
from pure RGI is still consistent with that using the RGI plus power correction ansatz. 
The overall shift of about 2% at Z energies indicates a possible influence of power 
terms. The spread of the unweighted mean as values is only of the order of 0.002. It 
indicates possible differences in power terms and uncalculated higher order contributions 
which are naively expected to be of the order ~ 0.002. The strongly differing second 
order contributions (compare the B/A values of the different observables in MS shown in 
Figure 14) as well as the differing NLLA expansions suggest differing higher order terms 
for the individual observables. 

As a further illustration of the validity of RGI one may try to calculate the size of the 
power terms observed in the MS analysis. The power term can be approximated as the 
difference of the RGI and the MS terms: 

= RRGl{Ecm) — Rusi^cm) (31) 
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ocsiMz) 




/I T'X 

(1 -T) 

(Major) 


u.iiyy±u.(juuz±(j.uuu2 

0.1166±0.0001±0.0003 


31.9/4(J 
10.1/14 


(U-raramctcr) 


U.izU5±U.UUUz±U.UUU2 


24.0/22 


{K /^vis) 
//\/r2 / Tp2 \ 

(M.^E/i5v\) 


0.1274±0.0005±0.0017 

U. 1zoz±U.U(J(Jd±(J.(JU(J5 

n 1 1 Qn_i_n nnn/ij_n nnnp; 
U.lloy±U.UUU4±U.UUU0 

0.1216±0.0004±0.0004 


46.9/26 

'7/^/1/1 

7.6/14 

n /I /I /I 
y.4/ 14 

8.1/14 


/ R \ 
X-'-'max/ 

(-Ssum) 


n lOQfi-un nnno-i-n nnn/i 
U. izou±U.UUUz±U.UUU4 

0.1194±0.0001±0.0004 


oz.o / zz 
21.5/22 


(EECyoo-llO") 

(EECsoo-isQo) 


0.1169±0.0002±0.0001 
0.1199±0.0001±0.0001 


18.8/14 
44.0/14 


(JCEFiino-i6o°) 


0.1155±0.0002±0.0001 


31.5/14 


weighted mean 
unweighted mean 


0.1195±0.0002±0.0010 
0.1201±0.0020±0.0003 





Table 9: Results for a pure RGI fit to a large set of measurements of different experi- 
ments [41]. For Ecra > Mz only Delphi measurements are included in the fit. The first 
error is the statistical uncertainty from the fit, the second one is the systematic uncer- 
tainty. For the mean values only the E-definition Jet Masses have been used and both 
EEC and JCEF have been omitted. For the definition of the mean values see section 2.4. 

for a fixed Ajy^ = 250 MeV. This value corresponds to the average as as obtained from 
the RGI analysis. As this ansatz is slightly energy dependent, Mz has been chosen as the 
reference energy. Figure 14 shows the comparison between the measured simple power 
correction terms and the calculated power term using pure RGI. The strong correlation 
between Ci/A and B/A for fully inclusive means (left side) is well represented by the 
above ansatz. It also holds reasonably well for the JCEF and the EEC. Here the ratio 
of the first and second order coefficient B/A can be adjusted by an appropriate choice of 
the X integral. 

The agreement of the measured values Ci/A and the RGI curve depends on the choice 
of Ajy^. A higher value of 280 McV as indicated by the upper edge of the grey band 
shown in Figure 14 requires the data points to move down about by the same amount, 
leading to a difference between the data and the curve of about the width of the grey 
band. 

Similarly to the simple power term Ci, the values of the Dokshitzcr- Webber model 
can be calculated and obviously agree reasonably with the fitted values, as shown in 
Figure 16. The measured ao values, show a scatter around a universal value anticipated 
by the Dokshitzer Webber model. Even this scatter is described reasonably by the RGI 
prediction. 

Finally a remark on the connection of RGI to optimised scales is in order here. Though 
reference [5] stresses the difference in philosophy of RGI perturbation theory and scale 
optimisation there is an obvious connection. 

RGI perturbation theory describes the data on fully inclusive shape means well with- 
out any renormalisation scale dependence whatsoever. This implies that a fit of the 
rcnormalisation scale using the standard 0{a^) expansion leads to a scale value where 
the scale dependence of the 0{a1) prediction vanishes. This is the PMS scale which is 



32 




^.7 
> 

o 

^0.6 
U 



0.5 



0.2 



0.1 



-0.1 



DELPHI 



• EEC 90"-160" 
A EEC 70"-110° 
+ EEC 30"-150" 

♦ EEC 40"-140" 
□ JCEF 110"-160° 

0.4 JCEF 90°-160° 



0.3 - 



=^ 



10 



15 



20 



B/A 



B/A 



Comparison of the results for Ci/A obtained from a fit with fixed 



Figure 14: 

ation of A^^^ between 0.220 GeVand 0.280 GeV. Left: fully inclusive observables. Right: 
different integration intervals of EEC and JCEF. Note that the data points are correlated. 



^MS ~ 0.250 GeV and the prediction from pure RGI. The shaded band shows the vari- 



close to the ECH scale. The above contiguity has been observed numerically [31] and 
found to be valid within errors. 



6.4 Measuring the QCD /^-function based on the RGI analysis 

Although RGI perturbation theory is based on the RGE the /3-function itself is not 
fixed by the RGI approach. As the /3-function for an observable R and the QCD coupling 
is identical in second order this implies that the /5-function, or rather the coefficients /5o 
and pi, can be directly inferred from the energy evolution of the event shape observable 
means (see Equation 27). As pi induces only a small correction ((9(4%)) and cannot be 
simultaneously determined with Pq given the current precision of the data, it has been 
fixed to the QCD expectation. Assuming QCD, the number of active flavours, Uf, can be 
determined as an alternative to determining Pq. In these cases /3o and pi were calculated 
from the QCD relations using the fitted nf and assuming A^c = Ca = 3. 

The approach allows /3o to be measured under the assumption that QCD is valid. In 
order to compare the result with predictions from e.g. super symmetric extensions of the 
standard model one of course cannot rely on the QCD calculation for f3i. Instead one 
needs to fit the full energy dependence. This analysis is described in section 6.5. 

The implicit Equation 27 can be fitted to data on shape observable means leaving A/j 
and Po (or nf) as free parameters. Note that this fit is independent of an Uf dependence 
entering implicitly through the second order perturbative coefficient B. This holds as the 
offset term /3o In Aji is unrestricted because Ar is a fit parameter and implies that a real 
measurement of /3o or nf can be performed. 

The results of the fits to Delphi data are given in Table 10. The systematic uncertain- 
ties of the data are fully propagated. The small uncertainty on the approximately inverse 
power behaved b-mass correction (see Figure 3) is given separately. Further systematic 
uncertainties due to the possible presence of small power corrections are not considered 
here. Contrary to the case of an measurement, this seems to be justified since the 
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Observable 


/3o 




X /ndf 


/ml j_ \ 

(Thrust) 


7.7±1.1±0.2±0.1 


4.6±1.3±0.3±0.1 


9.2/13 


(C — parameter) 

\^h/^vis / E det. 
{M!/El,, ) E def. 


7.8±l.U±U.3±U.l 


4.7±1.2±U.4±U.l 


7.2/13 


7 Kj-i K_Ln O-Un n 


A S-Li O-Un Q-un n 
4.o±i.y±u.o±u.u 


o.o/ io 


7.5±1.1±0.2±0.0 


5.0±1.3±0.3±0.1 


7.1/13 


(-^max) 


7.7±1.4±0.1±0.1 


4.7±1.9±0.1±0.2 


6.3/13 


(-Ssum) 


7.7±0.9±0.1±0.1 


4.8±1.2±0.1±0.2 


5.9/13 


(Major) 


8.0±1.1±0.1±0.1 


4.3±1.5±0.2±0.2 


9.2/13 


Weighted Mean 


7.7±0.9±0.1±0.1 


4.7±1.2±0.1±0.1 




Theory 


7.67 


5 





Table 10: Results for the fits of f3o or Uf for different observables with Delphi data. The 

first uncertainty is statistical, the second systematic, and the third indicates the uncer- 
tainty due to the b-mass correction. Here the minimal uncertainties of the individual 
observables are taken as uncertainties of the mean value. 

energy dependence is even less affected by the power correction. For all observables 
the results for or Uf, are fully consistent with the expected values from QCD. This 
observation further strengthens the confidence in RGl perturbation theory. 

In order to reduce the uncertainty on (3o or Uf more data, especially from low energy 
experiments, needs to be included to increase the fit range. As {M^/El^ receives large 
power corrections due to mass terms, in practice only (1 — T) is left over as a possible 
observable. The data given in [18,41,49] were used in a common fit. As systematic 
uncertainties of the low energy data cannot be controlled in detail for this fit, systematic 
and statistical uncertainties were combined leading to small y^/ndf values. Thus the 
uncertainties on the fit values already include systematic components. The fit describes 
the data very well in the energy range a/s = 14 to 200 GeV. The slope of the evolution 
is expected to change at the b production threshold at ^/s ~ 14 GeV. The central result 
of the fit is shown in Figure 17 and Table 11. This table also shows the results when the 
(1 — T) data of the experiments are fitted separately. All results are consistent with the 
QCD expectation within their fit error. In order to estimate a systematic uncertainty of 
/5o or rif , the measurements of one low energy experiment at a time have been left out from 
the combined fit. This leads to a maximum deviation from the central fit of A/9o = +0.11, 
Arif = -0.27 (excluding PLUTO) and A/3o = -0.08, Anf = +0.17 (excluding JADE). 
Adding the fit error and the systematic uncertainty of the Delphi data in quadrature to 
the above deviations from the central result leads to the final result: 

(3q = 7.86 ± 0.32 , 

Uf = 4.75 ± 0.44 . 

These have to be compared with the QCD expectation /3o = 7.67 or nf = 5. 

As a final cross-check of the RGI method one may try to determine the NNLO co- 
efficient p2 (see Equation 25) from a fit to data. Small values of p2 indicate a good 
approximation of the /3-function by the two loop expansion. 

As p2 is mainly sensitive to low energy data again only (1 — T) is used. Fixing h and 
Pi to the values expected from QCD a two parameter fit yields: 



as = 0.119 + 0.001 + 0.001 
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/)2 = 4±3±3.5 . 

The result for p2 is compatible with zero within its uncertainty, supporting the NLO 
approximation. 

Prom p2 the expected NNLO coefficient of the MS expansion can be deduced: 

C = 915 ± 40 . 

This indicates a sizable correction to the NLO prediction in the MS scheme. Even 
for p2 = the NNLO coefficient is large: C — 882. Using the Pade estimate for the 
parameter C leads to a value of p2 = —10.8 and a corresponding change of the strong 
coupling Aas — +0.003 with respect to the case when the NNLO term is omitted. 



6.5 Implications for light gluinos 

Knowledge of the f3 function implies model independent limits on the presence of 
hypothetical particles such as gluinos. Supersymmetric extensions of the standard model 
predict a weaker energy dependence due to the effect of gluino loops. In testing this 
hypothesis one should evidently not presume the QCD prediction for Pi as has been done 
for the Po measurement in the previous section. Instead the full logarithmic energy slope 

dlnQ2 i?2 ^"^ ) 

has to be measured. Evaluating this equation in second order implies a small dependence 
on R{Q) or the energy respectively. For our measurements the choice of the Z-mass as 
central scale is reasonable. For R — {1 — T)/A one gets: 

dR-l _Q0O 

dW " 

as the QCD prediction. Supersymmetry alters this prediction. The actual value of the 
logarithmic slope, however, depends on the mass of the additional particles. Gluinos in 

the mass range from 5 to 190 GeV arc excluded from direct searches [1]. Gluinos in the 
open mass range below 5 GeV can be considered massless compared to the energy range 
of this experiment. Therefore one expects (3q^^^ = '^'^ '^^^ and /Sf^*^^ = £L_1£^_ This 
yields the supersymmetry prediction for the logarithmic energy slope of: 

dW^'-'' ■ 



Experiment 


/5o 


x'/ndf 


nf 


Delphi 

L3 
JADE 
TASSO 
PLUTO 


7.7±1.1 ±0.1 

10.3±2.3 ±0.0 
7.8±0.6 ±0.2 
7.7±1.2 ±0.2 
8.3±1.3 ±0.2 


9.22/13 

0.34/4 

0.92/2 

0.02/2 

1.00/4 


4.6±1.3 ±0.1 
2.3±2.9 ±0.0 
4.8±0.7 ±0.2 
5.0±1.4 ±0.2 
4.2±1.5 ±0.2 


combined 


7.86±0.16±0.1 


25.7/37 


4.75±0.18±0.1 


Theory 


7.67 




5 



Table 11: Results for the fit of f3o for Thrust for different experiments, and as a combined 
fit. The combined fit includes additionally datapoints from AMY, TOPAZ, MK II and 
HRS. The first error is statistical, the second is the uncertainty of the b-mass correction 
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Measuring the slope by fitting the function l/{a + blnQ) to the data yields: 

JT) — 1 

^^8.70 ±0.35 . 

The error contains the statistical and systematic uncertainties in a similar way as dis- 
cussed for the (3q measurement in the previous section. The existence of light gluinos can 
be firmly excluded by our result. 



7 Summary and conclusions 

Event shape distributions and their means or integrals for the obscrvablcs 1 — T, 
C-Parameter, -Bmax, -Bsum, EEC, JCEF, the classical and the E- or p-scheme definitions 
of M^/Ey^g, Mg / E^^, have been presented in the wide range of hadronic centre-of-mass 
energies from 45 GeV to 202 GeV. Measurements at energies below ^/s — Mz were ob- 
tained from events containing a hard radiated photon. 

The data can be successfully described over the whole energy range using Monte 
Carlo fragmentation models or by analytical power correction models. The successful 
description of the EEC by a power correction ansatz is reported for the first time. 
Fitting the Dokshitzer-Webber power corrections to event shape distributions yields 
as = 0.1110 ± 0.0055rms ^-nd ao = 0.559 ± 0.073rms- The application to the mean 
values of event shapes yields as = 0.1217±0.0046rms and ao = 0.431 ±0.048rms respec- 
tively. A measurement of the shift of the Sudakov shoulder of the C-Parameter indicates 
an approximately constant power shift over the whole three jet range. 

In addition to the comparison of power correction models, the data for inclusive shape 
means of seven observables depending only on a single energy scale and integrals of the 
EEC and JCEF have also been compared to results of the so-called Renormalisation 
Group Invariant, RGI, perturbation theory with and without additional power terms. 
This method allows for a measurement of which is independent of any renormalisation 
scheme or scale. With respect to as the RGI prediction is equivalent to the ECH scheme. 

It has been observed that RGI perturbation theory is able to describe the energy evo- 
lution of these data consistently with a single value of the strong coupling parameter: 
as{Mz) — 0.1201 ± 0.0020j?MS- Within RGI perturbation theory there is no need for 
power corrections. The small R.M.S. of the as values obtained from the different observ- 
ables indicates an improved decription of the data by RGI perturbation theory compared 
to the standard MS treatment. Furthermore it serves as an important consistency check 
of the method. 

The most important single result of the analysis is the measurement of the /3-function 
of strong interactions. For the leading coefficient , assuming QCD, the average of seven 
observables as measured by DELPHI is: 

/?o = 7.7 ±0.9 ±0.1 , 

where the first uncertainty is statistical and the second systematic. This corresponds to 
the number of active flavours: 

71/ = 4.7 ± 1.2 ±0.1 . 

The systematic uncertainty accounts for experimental uncertainties as well as for the 
uncertainty induced due to the correction for b hadron decays. Inclusion of additional 
low energy data for the observable (1 — T) yields the result: 

Po = 7.86 ± 0.32 , 
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Uf = 4.75 ± 0.44 . 

Here the uncertainty includes the sources mentioned above as well as an estimate of the 
systematic error induced by the inclusion of the low energy data. 

Within RGI this quantity can be derived without any renormalisation scheme depen- 
dence. Power corrections have been found to be negligible. The precision of this result 
is greatly increased compared to previous determinations of /Sq from event shape observ- 
ables [49] as well as to a determination based on the most precise measurements of as[50]: 
Po — 7.76 ± 0.44. Further reduction of the uncertainty is to be expected from a proper 
combination of the results of the LEP experiments. The analysis should then also be 
repeated for observables other than (1 — T). 

Fitting directly the logarithmic energy slope yields 

JTD— 1 

j^ = W = 8.70 ±0.35 . 

This measurement excludes light gluinos in the open mass range below 5 GeV [1] in a 
model independent way. 

Provided the possible presence of power terms in the event shape means is clarified 

by future studies, the good stability of the results for for a larger number of event 
shape observable means observed when using RGI perturbation theory may indicate the 
possibility of improved measurements of the strong coupling as- 
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